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ABSTRACT 


Simple  acoustic  source  theory  has  been  applied  to  determine 
the  maximum  overpressure  obtainable  by  the  deflagration  of  non- 
spherical  clouds.  In  three  dimensions  overpressure  is  generated 
not  by  the  rate  of  energy  addition  but  by  the  first  time  deriva¬ 
tive  of  the  rate.  Because  of  this,  deflagrative  combustion  of 
edge-ignited  clouds  produces  markedly  less  overpressure  than  cen¬ 
tral,  spherical  ignition.  Examples  are  presented  for  three  non- 
spharical  cloud  and  igniter  geometries.  The  implication  is  that 
even  high  velocity,  subsonic  combustion  waves  cannot  produce 
damaging  blast  waves  and  that  some  type  of  supersonic  combustion 
or  massive  flame  acceleration  is  required  if  a  damaging  blast  wave 
is  to  be  produced. 
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INTRODUCTION 


The  determination  of  the  blast  wave  produced  in  free  space 
by  the  nonspherical  combustion  of  a  nonspherical  vapor  cloud  has  not 
been  examined  in  any  detail ,  either  analytically  or  numerically. 

The  problem  is  important  because  at  the  present  time  spherical  theory 
is  being  used  to  estimate  the  overpressure  for  clouds  of  arbitrary 
shape  and  no  one  knows  how  conservative  this  approach  is.  This  paper 
applies  simple  acoustic  source  theory  to  the  deflagrative  combustion 
of  a  few  simple  free  cloud  shapes  to  determine  a  first  cut  estimate 
of  the  maximum  overpressure  that  one  could  expect  as  a  function  of 
cloud  size  and  shape. 


THE  SIMPLE  ACOUSTIC  SOURCE 

Stokes^  first  showed  that  a  simple  source  of  mass  at  a 
point  m(t)  (Units:  Kg/ sec)  generates  a  sound  wave  in  three  dimen¬ 
sions  which  has  an  overpressure  which  is  proportional  to  m(t). 
Specifically,  as  Lighthill^  states, 

P  -  PQ  =  fijt  -  ~J/4*r 

where  t  -  ~~  replaces  t  because  the  wave  is  propagating  away  from 
*o 

the  source  region  at  the  velocity  of  sound,  aQ.  Lighthill  also 
states  that  at  some  distance  from  the  source  the  sound  will  appear 
to  emanate  from  a  point  (i.e.,  the  wave  will  be  very  close  to  spheri¬ 
cal)  if  the  quantity  (ur/aQ)a  '«  1.  Here,  w  is  the  circular 
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« 


frequency  of  the  source  of  sound  and  r  is  a  characteristic  radius 
of  the  source  region.  Finally,  for  a  source  which  has  a  finite 
duration  /m(t)dt  «  0,  and  therefore  /(P-PQ)dt  =  0  in  the  wave. 

This  means  that  the  acoustic  pulse  from  a  source  of  finite  dura¬ 
tion  must  have  a  negative  phase  impulse  equal  to  the  positive 
phase  impulse.  This  should  be  contrasted  to  the  pulse  emitted  by 
a  finite  duration  source  in  a  strictly  one-dimensional  channel. 

In  this  case  the  acoustic  overpressure  is  proportional  to  m(t)  and 
the  overpressure  in  the  pulse  is  always  positive  everywhere.  ^ 

Now  consider  the  deflagrative  combustion  of  an  unconfined 
cloud  of  arbitrary  shape.  For  simple  source  acoustic  behavior  the 
effective  rate  of  mass  addition  m(t)  may  be  replaced  by  an  effec¬ 
tive  rate  of  volume  addition  multiplied  by  the  initial  gas  density, 
that  is,  m(t)  =  PQv(t}. 

However,  for  deflagrative  combustion 


m(t) 


Po^C*)  °  P0 


5F 


Af(t)) 


(2) 


Where  Su(t)  is  the  effective  normal  burning  velocity,  A^(t)  is  the 

effective  frontal  area  of  the  flame  and  the  ratio  (V. -V  )/V  is  the 

new  volume  produced  per  unit  volume  of  gas  burned  at  constant  pres- 
* 

sure.  Using  the  notation  of  Strehlow  et_  al .  ^  (V^-V^/V^  ■  §/yd. 
Here,  $  is  an  effective  dimensionless  heat  addition.  Its  relation 
to  the  q  of  Strehlow  et  al.  is  given  in  Appendix  A.  Therefore, 
using  the  definition  of  the  velocity  of  sound  in  an  ideal  gas  and 
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combining  Eqs.  (1)  and  (2)  yields 

•  V1” 

where  P  is  the  acoustic  overpressure.  Equation  (3)  is  a  signifi¬ 
cant  relationship.  It  states  that  the  acoustic  overpressure  gen¬ 
erated  by  deflagrative  combustion  in  three  dimensions  is  given  by 
the  time  rate  of  change  of  the  rate  of  energy  addition  rather  than 
the  rate  of  energy  addition  as  it  is  in  strictly  one-D  flow. 

The  deflagrative  combustion  of  a  cloud  generates  a  sound 
wave  of  very  low  frequency.  The  period  of  the  sound  source  in 
this  case  is  given  by  x  =  2ir/«  where  T  is  the  burn  time.  Thus  the 
requirement  that  the  source  region  is  acoustically  simple  becomes 

l27TCrB/T)/a0J2  «  1 

where  r  is  a  major  dimension  of  the  cloud  and  approximates  the 

ID 

distance  traversed  by  the  deflagration  wave.  Even  with  distortion 
and  motion  due  to  flame  propagation  r/x  “  Su  and  the  simple  source 
requirement  becomes  4ir2Mgu  «  1.  If  the  quantity  is  to  be  less 
than  0.1  for  no  blast  wave  distortion  Su  will  be  about  17.5  m/sec 
before  deflagrative  behavior  produces  significant  distortion  of 
the  blast  wave.  Even  for  significantly  higher  values  of  Su  the 
first  effect  will  be  only  to  distort  the  shape  of  the  blast  wave. 
Higher  overpressures  than  estimated  here  will  not  occur  until  Su 
becomes  quite  large.  This  limitation  of  the  theory  has  not  been 
quantified  as  yet. 
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APPLICATION  IN  SPHERICAL  GEOMETRY 
Two  classical  problems  have  been  treated  in  spherical  ge¬ 
ometry.  These  are  the  central ly-ignited,  constant -velocity 
flame and  the  Taylor  piston  problem.  ^  For  a  constant-velo¬ 
city  flame  Eq.  (3)  reduces  to  the  relation 

p.  Jfu  ^ 

4ira*r  dt 

where  r  is  the  distance  from  the  source  to  the  observer.  For  a 
growing  spherical  flame  ball 

A^  =  4irr* 


and 


drf  Vh 

dT8SssS«r 


(4) 


Substituting  yields 


2$S2r. 

P  -  AhJ. 

.  u 


where  the  r  in  the  denominator  must  be  equal  to  or  greater  than  r^. 
If  we  let  it  be  equal  to  r^,  and  if  we  substitute  the  definition  of 
Q  given  above,  we  obtain  an  expression  for  the  acoustic  pressure  on 


the  unbumed  side  of  the  flame 


2Y 


vl 

fvj 

1  -  — 
v. 

b 

V 

w 

Mi 


Su 


(6) 


which  is  identical  to  Eq.  (19)  of  Strehlow  et  ad. 


(S) 


If  we  write  =  Sgt  for  the  flame  propagation  in  Eq.  (5) 
and  recall  that  away  from  the  flame  t  becomes  jt  -  Eq.  (5)  may 
be  written 


P  =  2y  1 


V  fv. ) 3  fant 
.  Ji  _k  M3  _5L  _  i 
V.V  Su  r 

bj  uj 


However,  the  transformation  from  a  flame  to  a  Taylor  spheri¬ 


cal  piston  is 


V  }J/*  V 

M  *  1  -  ~  tr-  M_ 

p  Vfe  Vb  Su 


Therefore,  Eq.  (7)  reduces  to 


P  *  2y 


fa  t  } 

M3  - - lj 

°  Pi  r  J 


where  the  range  on  r  is  S^t  <  r  <  aQt.  Equation  (8)  is  identical 
to  the  Taylor  solution  for  a  spherical  piston  as  0  (Eq.  (13), 

Strehlow  et^  <al.  ^ ) . 

We  note  that  in  these  equations  the  rate  of  flame  area  growth 
is  proportional  to  r  and  the  rate  of  decay  of  a  spherical  acoustic 
pulse  is  proportional  to  r"J.  Thus,  for  the  spherical  case  the 
acoustic  pressure  remains  constant  at  the  flame  front  and  the  wave 
that  is  generated  during  flame  motion  is  self-similar  in  r/t.  This 
is  a  fundamental  acoustic  property  of  a  constant -velocity  flame  that 
is  completely  surrounded  by  combustible  material  when  it  burns.  It 
has  been  verified  by  numerous  theoretical  and  numerical  investiga¬ 
tions.  (3'4) 


APPLICATION  TO  FREE  CLOUDS 

Ke  restrict  ourselves  to  clouds  of  arbitrary  shape  in  free 
space.  To  simplify  the  treatment,  we  will  consider  free  clouds 
which  have  a  horizontal  plane  of  symmetry  for  cloud  shape,  cloud 
distortion  during  deflagration  and  deflagration  wave  shape.  Since 
we  now  wish  to  simply  examine  the  implications  of  Eq.  (3)  for 
general  geometries,  we  for  the  moment  restrict  ourselves  to  con¬ 
sidering  only  constant  velocity  flames.  This  means  that  Eq.  (3) 
reduces  to  Eq.  (4)  where  is  an  effective  normal  burning  velocity 
and  r  is  the  distance  of  an  observer  from  the  cloud  center. 

In  general  the  evaluation  of  the  dA^/dt  term  in  Eq.  (4)  is 
complicated  because  during  deflagrative  combustion  the  cloud  ahead 
of  the  flame  will  move  and  distort  with  time.  This  distortion  will 
cause  A^  to  be  a  different  function  of  time  than  that  calculated 
from  purely  geometric  considerations  using  a  constant  velocity  flame 
propagating  through  a  stationary  cloud. 

In  order  to  determine  how  to  best  approximate  the  effect  of 
such  distortion  for  arbitrary  ignition  of  an  arbitrary  cloud  one 
must  first  examine  flow  behavior  in  two  simple  limit  geometries; 
strictly  one-dimensional  flow  caused  by  deflagration  and  volumetric 
heat  addition  in  three  dimensions. 

The  wave  diagram  for  a  strictly  one-D  flow  driven  by  a  con¬ 
stant-velocity  deflagration  wave  is  shown  in  Fig.  1.  If  we  model 
the  combustion  of  a  stoichiometric  methane-air  mixture  as  proposed 
in  Strehlow  et  al.,^  the  gas  to  the  right  of  the  flame  and  to  the 
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left  of  the  contact  surface  (regions  0,  1  5  3)  can  be  assumed  to 

have  a  heat  capacity  ratio  yo  -  Yj  =  "  1.4,  while  the  gas  in 

region  2  has  a  =  1.202.  Additionally,  the  effective  value  of 

dimensionless  energy  addition  is  q  =  8.430.  This  means  (from 

Appendix  A)  that  q  =  9.242  and  that  V^/Vu  =  7.602  for  constant 

pressure  combustion.  Furthermore,  for  this  system  the  Chapman- 

Jouguet  deflagration  Mach  number  is  0.16631. 

The  properties  of  the  wave  system  of  Fig.  1  were  calculated 

using  the  shock  relationships  from  Liepman  and  Kosliko^  for  the 

shocks  Sj  and  S^,  the  deflagration  relationship  from  Strehlow  et 
(31 

al. v  J  for  the  flame  and  the  contact  surface  requirement  that 

U  =  U  and  P_  =  P_.  The  results  of  this  calculation  that  are 
p2  p3  2  3 

important  to  this  paper  are  summarized  in  Fig.  2. 

Figure  2  shows  that  particle  velocities  generated  by  flame 
propagation  are  relatively  high  and  that  the  apparent  flame  speed 
relative  to  an  observer  is  from  4.3  to  9.3  times  the  actual  normal 
burning  velocity.  Also  it  is  less  than  the  space  velocity  of  a 
spherical  flame  if  the  one-dimensional  flame  is  traveling  at  less 
than  about  0.9  of  its  maximum  (CJ)  velocity.  Furthermore,  the  ratio 
| Up  /Up  |  lies  between  1.0  and  1.4  in  this  same  velocity  range. 

This  last  observation  simply  means  that  even  relatively  high 
speed  flame  propagation  from  a  free  surface  in  one-dimensional 
space  displaces  the  gas  rather  uniformly  to  the  front  and  back. 

This  is  to  be  contrasted  to  the  behavior  of  a  one-dimensional  CJ 
detonation,  which  displaces  the  gas  primarily  in  1  he  direction  of 
detonation  propagation. 
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We  now  wish  to  determine  what  happens  if  we  relax  the  con¬ 
finement  normal  to  the  flame  propagation  direction  for  the  case  of 
deflagrative  combustion.  To  estimate  this  effect  we  assume  that  if 
energy  is  added  relatively  slowly  to  a  gaseous  volume,  the  volume 
will  expand  an  equal  distance  in  all  directions.  During  flame 
propagation  through  a  cloud  the  flame  can  be  assumed  to  be  a 
relatively  thin  sheet  of  energy  addition  when  compared  to  the  cloud 
dimensions.  In  other  words  the  flame  can  be  assumed  to  occupy  a 
volume  AAL  when  AL  «  /k.  Under  these  circumstances  one  would 
expect  the  major  expansion  to  occur  in  a  direction  normal  to  the 
local  orientation  of  the  flame  even  though  the  pressure  is  propa¬ 
gating  almost  equally  in  all  three  dimensions. 

This  assumption  has  beer,  qualitatively  verified  experimen¬ 
tally  by  observing  flame  propagation  in  a  propane-air  elongated 
pancake  cloud  above  a  flat  plate  of  about  50  mm  height  and  0.6  m 
x  0.3  m,  edge  ignited  at  the  center  of  the  0.3  m  edge. ^  The 
flame’s  apparent  propagation  speed  was  about  3.5  times  the  normal 
burning  velocity,  and  its  luminous  height  was  about  equal  to  the 
initial  cloud  height.  Furthermore  the  flame  flashed  to  encompass 
an  area  considerably  larger  than  the  original  lateral  extent  of  the 
cloud.  The  assumption  has  also  been  reasonably  well  verified  by 
numerical  calculations  made  by  Dr.  Lon  Hazclmanv  ‘  at  Lawrence 
Livermore  Laboratories  using  a  2  D  hydrocode.  He  calculated  flow 
and  flame  propagation  in  an  open  2  D  channel.  His  calculations 
yielded  a  flame  height  almost  twice  the  initial  cloud  height  and  an 
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apparent  flame  speed  about  three  times  the  normal  burning  velo¬ 
city.  However,  his  numerical  flame  had  a  thickness  only  one-third 
of  the  cloud  height  and  this  should  lead  to  more  vertical  dis¬ 
tortion. 

These  results  indicate  that  for  any  arbitrary  free  cloud  and 
igniter  geometry  in  which  the  flame  is  not  completely  surrounded  by 
a  combustible  mixture  one  may  make  the  following  conservative  as¬ 
sumptions  about  flame-induced  distortion. 

1.  The  burning  velocity  that  should  be  used  is  S^.  This  is 
because,  irrespective  of  the  motion  of  the  flame  induced  by  the  flow, 
the  entire  cloud  will  burn  in  a  time,  T,  given  by  S^r  =  rm  where 

is  the  distance  from  the  igniter  to  the  farthest  point  in  the  ini¬ 
tially  quiescent  cloud. 

2.  The  increase  of  flame  area  in  a  free  cloud  due  to  dis¬ 
tortion  of  the  cloud  by  flame-induced  flow  will  never  exceed  the 

value  V,  /V  calculated  for  isoboric  combustion.  This  is  because 

b  u 

this  is  the  value  of  distortion  of  the  unburned  material  produced 
during  combustion  of  a  completely  surrounded  (spherical)  flame,  and 
relief  due  to  product  gas  motion  away  from  the  flame  can  only  re¬ 
duce  this  factor. 

Thus  the  general  procedure  for  determining  a  first  conserva¬ 
tive  approximation  to  tho  effect  of  cloud  shape  and  igniter  loca¬ 
tion  on  the  overpressu; e-time  behavior  in  the  blast  wave  when  the 
igniter  is  not  completely  surrounded  by  the  cloud  is  as  follows: 

1.  Assume  a  value  of  Su> 
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2.  Assume  a  cloud  geometry  and  igniter  location. 

3.  Calculate  dA^dt  for  a  flame  traveling  at  velocity 
Su  through  the  quiescent  cloud. 

4.  Correct  for  distortion  by  multiplying  by  V^/V^ 


Examples  for  Specific  Geometries 

a.  Central  (line)  ignition  of  a  pancake  cloud. 
The  flame  area  is  given  by 


Af  =  2rrrfH 

where  r^  is  the  flame  radius  and  H  is  the  height  of  the  pancake. 
Therefore 


and 


dr 

2,111  a r =  2,THSu 


p 

max 


$s2h 

^  u 

2a*R  . 
o  obs 


(9) 


where  RQ^s  is  the  distance  of  an  observer  from  the  center  of  the 
cloud. 

b.  Edge  (line)  ignition  of  a  pancake  cloud. 

The  area  of  a  flame  propagating  from  an  edge  of  a  cloud 
of  radius  R  is  given  in  terms  of  the  distance  propagated  by  the 
formula 

r- 

Af  -  2Hrfcos_1  ^ 

Now 


Therefore 


or 


max 


$S2H 

u 

4a*R  . 
o  obs 


(10) 


which  is  exactly  half  that  predicted  for  a  centrally-ignited 
pancake. 


c.  End  (point)  ignition  of  an  ellipsoid  (cigar-shaped)  cloud 
Assume  a  flat  flame  of  location  r^  propagates  from  one 
end  of  a  long  cloud  with  a  minor  radius  of  rotation  d  and  major 
radius  b. 

Flame  area  is  therefore 


or 


2rrr 


d2 
f  b* 


Now  let 


•  b  •  ■  L/2  <  r  i  L/2  «  b  and  D  *  2d 


where  L  is  the  length  of  the  ellipsoid. 
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Therefore 


ind 


P 


-qS2D2  fv. 
u  b 

2a*rL2  rf  V 
° 


P 

max 


D 

L 


Cll) 


Now  assume  that  Su,  aQ  and  (V^/V^  are  independent  of 
cloud  shape  and  that  the  observer  is  at  a  fixed  distance  from  the 
center  of  a  cloud  of  fixed  volume  but  of  these  different  shapes. 
Call  the  volumes  *»  Vp  *  Vc,  where  the  subscripts  s,  p  and  c  re¬ 
fer  to  spherical ,  pancake  and  cigar  shaped  clouds,  respectively. 

Also  define  the  aspect  ratio  of  the  cloud  /4l  as  /H  =  2R  /H  and 

P  p  P 

=  Dc/Lc.  Using  this  notation  we  obtain 


(12) 


and 


c  r* 

Combining  Eqs.  (5),  (9),  (10)  and  (11)  with  Eqs.  (12)  and  (13) 


(13) 


yields 


(14) 


i 

i 

l 

i 

i 
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i 
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and 


IS 


max 


1»/» 


t  *  *)  * 


(15) 


where  the  factor  ♦  1  ■  V.  /V  arises  because  we  are  considering 

the  initial  cloud  diameter  and  Eq.  (5)  is  written  in  terms  of  the 

final  radius.  Eq.  (14)  is  for  central-line  ignition  of  a  pancake: 

for  edge-line  ignition  (P  /P  )  is  one  half  of  this  value. 

c  s  max 

Figure  3  is  a  plot  of  Eqs.  (13)  and  (14).  It  shows  that  the 
deflagration,  after  edge  ignition  of  a  large  aspect -ratio  cloud, 
produces  much  lower  overpressures  than  central  ignition  of  a  spheri¬ 
cal  cloud.  It  should  be  pointed  out  that  point  ignition  of  the 
edge  of  a  free  spherical  cloud  =  1  would  yield  a  higher  over¬ 
pressure  than  calculated  with  Eq.  (14)  because  the  derivation  of 
Eq.  (11)  assumed  a  flat  flame.  In  reality  one  would  expect  a 
flame  of  roughly  spherical  shape  to  propagate  away  from  a  point 
source  and  the  initial  rate  of  flame  area  increase  would  be  larger 
in  this  case  than  for  a  flat  flame. 


EFFECT  OF  CLOUD  SIZE 

Equation  (3)  can  be  used  to  give  an  order  of  magnitude  esti¬ 
mate  of  flame  areas  and  accelerations  that  are  necessary  for  a  def¬ 
lagration  to  produce  a  damaging  blast  wave.  If  we  define  the 
threshold  of  damage  to  be  0.1  Bar  and  assume  that  this  level  of 
overpressure  is  to  be  produced  100  m  from  the  cloud  center  in  an* 
atmosphere  with  a  velocity  of  sound  of  3S0  m/s,  Eq.  (S)  reduces  to 


1.6  x  10' 


A-  -rr—  ♦  S  t— 
f  dt  u  dt 


Table  I  was  constructed  from  Eq.  (16)  under  the  two  limit  assump¬ 
tions  that  (1)  dSu/dt  =  0  and  (2)  dAf/dt  =  0. 


TABLE  I 


(1)  dSu/dt  =  0  (2)  dAf/dt  «  0 


Su,  m/s 

dA, 

air*  ■*/* 

A£,  m* 

as 

dt** 

1 

1.7  x  107 

100 

1.7  x  10s 

10 

1.7  x  10* 

10,000 

1.7  x  10* 

100 

1.7  x  10s 

1,000,000 
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Notice  from  Table  1  that  even  for  very  high  flame  veloci¬ 
ties  the  rate  of  flame  area  increase  must  also  be  very  high  if  even 
a  weak  blast  wave  is  to  be  generated.  On  the  other  hand,  only 
extremely  large  initial  flame  areas,  exhibiting  very  large  flame 
acceleration  over  their  entire  frontal  area  are  necessary  to  pro¬ 
duce  a  weak  blast  wave  from  deflagrative  combustion.  Both  of  these 
observations  imply  that  extremely  large  clouds  are  required  if  one 
is  to  produce  significant  overpressure  by  deflagrative  combustion 


alone. 


DISCUSSION  AND  CONCLUSIONS 

The  application  of  simple  source  acoustic  theory  to  defla¬ 
grative  combustion  of  an  unconfined  cloud  shows  that: 

1)  It  is  very  difficult  to  produce  a  damaging  blast  wave 
by  deflagrative  combustion  once  the  flame  is  not  completely  sur¬ 
rounded  by  a  combustible  mixture. 

2)  The  maximum  overpressure  produced ,  other  circumstances 
being  equal,  is  proportional  to  the  ratio  of  the  minor  dimension 
of  the  cloud  to  the  distance  to  the  observer. 

3)  Blast  pressure  is  rather  uniformly  distributed  in  all 
directions,  i.e. ,  the  blast  is  roughly  spherical. 

4)  The  cloud  must  be  very  large  if  a  damaging  blast  wave 
is  to  be  produced. 

5)  Spherical  flame  propagation  calculations  such  as  those 
of  Kiuhl  et  al.  and  Strehlow  et  al.  greatly  overestimate  blast  pres¬ 
sures  from  deflagrative  combustion  following  edge  ignition  of 
clouds  with  large  aspect  ratios. 

Fran  various  accident  accounts  of  the  sequence  of  events 
that  led  to  the  production  of  a  damaging  blast  wave  after  delayed 
ignition  of  a  massive  spill  of  combustible  material  it  appears  that 

1.  There  is  a  threshold  spill  size  below  which  blast  damage 
does  not  occur.  Gugan's^  documentation  of  incidents  shows  that 
blast  damage  has  been  observed  for  spills  of  less  than  2000  Kg  but 
more  than  100  Kg  only  for  the  fuels  H_,  Hj-CO  mixture,  CH^  and 
CjH^.  Blast  damage  has  been  recorded  only  for  spills  greater  than 


16 


2000  Kg  for  all  other  fuels. 

2.  In  the  majority  of  cases  where  blast  damage  occurred, 
fire  was  present  for  a  considerable  period  before  the  blast  oc¬ 
curred. 

3.  In  many  cases  damage  is  highly  directional. 

These  observations  when  coupled  with  the  results  of  simple 
source  acoustic  theory  for  deflagrative  combustion  lead  to  the 
following  conclusions. 

1.  There  should  be  a  size  threshold  below  which  blast  will 
not  occur  as  long  as  ignition  is  "soft",  i.e.,  does  not  directly 
trigger  detonation. 

2.  The  fact  that  fire  is  present  early  after  ignition  in¬ 
dicates  that  massive  flame  accelerations  are  necessary  to  lead  to 
blast  wave  formation.  Since  the  flame  must  have  burned  through 
the  cloud  edge  by  the  time  the  blast  is  produced,  simple  acoustic 
source  theory  must  be  operative  for  even  high  deflagration  velo¬ 
cities.  Thus,  the  blast  must  arise  from  some  sort  of  effectively 
supersonic  combustion  process  or  from  very  rapid  increases  in  ef¬ 
fective  surface  area  of  the  flame. 

3.  Simple  acoustic  source  theory  for  deflagrative  combustion 
shows  that  deflagrative  processes  per  se  cannot  produce  highly 
directional  effects.  However,  it  is  well  known  that  detonative 
combustion  of  a  cloud  does  produce  highly  directional  blast  wave 
effects. 
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APPENDIX  A 

THE  GENERALIZED  HEAT  ADDITION  HUGONIOT 


Strehlow  e£  aK  ^  have  shown  that  a  generalized  heat  addition 
Hugoniot  quite  accurately  fits  the  real  Hugoniot  for  hydrocarbon- 
air  combustion  over  the  range  of  interest  for  vapor  cloud  explosions. 
Briefly  they  assume  that 

H  =  C  T 
o  p 
*o 

and 

h.  =  C  T  -  Q 
1  >>i 

where  Q  is  the  energy  added  at  the  deflagrative/detonative  discon¬ 
tinuity  and  C  and  C  are  different  heat  capacities  for  the  re- 
po  P1 

active  and  product  gases  respectively. 

Using  this  approach  they  define  a  dimensionless  energy  addition 


as 

q  -  CT,  -  DQ/P0V0 

Thus  the  volume  change  for  constant  pressure  combustion  is 
given  by  the  expression 


therefore 

/ 


iL 

yi 


V. -V 
b  u 


A* 

JL 

Yl 


v  y 

u  'o 


Now  define 
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for  ease  of  nondinensionlizing  the  acoustic  problem.  This  yields 


$ 


WV  Yo  - 

+  u-  q 


Vv1J  Yi 

For  stoichiometric  methane-air  mixtures^^ 


Yo  *  1.4 


Yj  =  1.202  and  q  = 


8.430 


Therefore 


Thus 

and 


=  £  =  6.6015 
u  To 

q  =  9.242 


»  7.6015 


for  a  stoichiometric  methane  cloud. 
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Figure  2.  Flow  velocities  and  apparent  speeds  associated 
with  the  strictly  one-dimensional  deflagration 
combustion  illustrated  in  Figure  1. 
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Figure  3.  Effect  of  the  aspect  ratio,  on  the  maximum 
blast  wave  pressure  rise  for  the  deflagrative 
combustion  of  pancake  and  cigar  shaped  clouds. 
Cloud  volume,  normal  burning  velocity  and  ob¬ 
server  distance  from  cloud  center  all  assumed 
to  be  constant  from  cloud  to  cloud. 


